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Dynamic quantities related to hysteresis have been measured in micro- 
magnetic simulations of single-domain nanoscale magnets at nonzero tempera- 
ture. The hysteresis-loop area and magnetization-field correlation display the 
characteristics of resonance, and the resonance frequency is found to be tem- 
perature dependent. The period-averaged magnetization displays symmetry 
breaking at high frequencies. 
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The common theme in hysteretic systems is a nonlinear, irreversible response that lags 
behind an applied force. While different physical mechanisms may cause this behavior, 
in many interesting situations it results from a system with two local free-energy minima 
separated by a free-energy maximum. For instance, the crystalline or shape anisotropy 
of a magnetic particle may create a barrier that significantly interferes with the reversal 
mechanism whereby the magnetization changes from the unfavorable alignment antiparallel 
to the applied field to the equilibrium alignment parallel to the field. Historically, hysteresis 
experiments have been conducted under quasistatic conditions that minimize thermal and 
dynamic effects. With new technology making applications of magnetization reversal of 
nanoscale magnets on the nanosecond time scale important, a thorough understanding of 
these effects is needed. Here we present results from dynamic micromagnetic modeling of 
hysteresis in nanoscale pillars at nonzero temperatures. 
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Our numerical model is based on iron nanopillars fabricated using STM-assisted chemical 
vapor deposition |]||. Following the one-dimensional model for magnetic pillars discussed 
by Boerner and Bertram ||, the nanoscale magnets are decomposed into a one-dimensional 
stack of 17 cubes with side 7.2 nm. Each cube has a magnetization density M(rj,t) and 
a local magnetic field H(r^t), with rj the center of the i-th cube. The damped preces- 
sional motion of each magnetization vector is given by the Landau-Lifshitz-Gilbert (LLG) 
equation (||| 

= -J2—M(n) x ( H( r< ) - — M( ri ) x MrA) , (1) 

where the electron gyromagnetic ratio is 7o=1.76xl0 7 Hz/Oe ||, and a=0.1 is a phenomeno- 
logical damping parameter. The numerical details have been given elsewhere Pa- 
rameters corresponding to bulk iron were used: saturation magnetization density M s = 
1700 emu/cm 3 and exchange length l e = 3.6 nm. The local magnetic field H(r», i) is a linear 
superposition of terms corresponding to exchange, magnetostatic interactions, the applied 
field, and a stochastic field representing thermal fluctuations. The latter term is governed by 
the fluctuation-dissipation theorem, for which we use the form derived for isolated magnetic 
particles |||§. The stochastic integration [g] of Eq. (JJ) was performed using a first-order 
Euler scheme with At=5xl0~ 14 s. 

Thermal fluctuations are essential in this model because the free-energy barrier can be 
surmounted using thermal energy momentarily "borrowed" from the surroundings. For the 
spatially nonuniform reversal mechanisms considered here , this barrier crossing occurs 

at the ends of the pillar. The observed hysteretic behavior depends on the interplay of three 
time scales: the mean nucleation time, the time required for the subsequent growth of the 
nucleated volumes to fill the particle, and the period of the applied field. 

To simulate hysteresis loops, an external field H = zH cos (2irvt) is applied, with 
i/o=2000 Oe. The results are presented for the reduced quantities m(i), the z-component of 
the magnetization for the entire magnet, and h(t), the applied magnetic field, both normal- 
ized to have maximum values of unity. 
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The correlation between m and h is measured by 
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m(t)h(t)dt , 



(2) 



where the integral is taken over one hysteresis loop. The results are presented vs. v in 
Fig. 1 for T=20K and 100 K. The bars shown for T=100K are not error bars, instead they 
represent the standard deviation of the probability density for B. An example of m(t) is 
shown as the solid curve in the inset of Fig. 1, along with h(t)/ 10 which appears as the dotted 
curve. In the high-frequency limit the magnetization usually does not switch, and B is near 
zero. At low frequencies the magnetization switches nearly every period, and B approaches 
unity. The quantity B corresponds to the reactive part of the nonlinear response function || , 
and the lowest-frequency zero crossing of B in Fig. 1 is associated with resonance. 
The dissipative part of the nonlinear response is the hysteresis-loop area, 



The measured average values of A are shown vs. v in Fig. 2 for both temperatures. At high 
frequencies the magnetization is unable to respond to the applied field, and the hysteresis 
loop area is nearly zero as most of the time m(t) oscillates with a small amplitude about one 
or the other of the degenerate values of the zero-field magnetization. In the low-frequency 
(or quasi-static) limit the magnetization usually switches at small values of h(t), and as a 
result A is again near zero. At intermediate frequencies A displays a maximum near the 
resonance frequency. The maximum in A and the low-frequency zero crossing in B agree 
to about 15%. The temperature dependence of the results is mostly due to the different 
average nucleation times associated with the different temperatures. Similar simulations in 
two-dimensional Ising magnets have shown that the phenomena seen here are associated 
with stochastic resonance 0. 

Another interesting quantity to consider is the period-averaged magnetization, 
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which corresponds to a dynamic order parameter for the system. Histograms are shown for 
several applied-field frequencies at T=100K in Fig. 3. At high frequencies the histogram 
has two maxima, corresponding to the magnetization being oriented predominantly in the 
positive or negative direction during the entire hysteresis loop. At low frequencies the 
histogram is centered around zero since the magnetization switches during every half-period 
of the applied field. Previously, a phase transition in this dynamic order parameter has 
been seen in Ising models of nanoscale mag nets when many nucleation events 

contribute to each reversal of the magnetization. In these models, this phase transition 
occurs near the frequency where B(v) has its high-frequency zero crossing. However, the 
detailed analysis and large statistical sampling required to make such a determination has 
not been performed for the present simulations. 

In summary, micromagnetic simulations of nanoscale magnets at nonzero temperatures 
and gigahertz frequencies display the characteristics of temperature-dependent resonance. 
Future work exploring the connections to stochastic resonance and dynamic phase transitions 
is planned. 
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FIG. 1. The correlation B, defined in Eq. (|2|), between the normalized magnetization m and 
the normalized applied field h, shown vs. the applied field frequency v, for T=20K and 100 K. 
The bars are not error bars; instead, they represent the standard deviation of B at T=100K. The 
standard error in the measured mean is smaller than the symbol size. As B is the reactive part of 
the nonlinear system response, the first zero crossing is taken as a sign of resonance. A particular 
example of m(t) (solid curve) and h(t)/W (dotted curve) for T=100K and v = 0.2 GHz is shown 
for several hysteresisdoop periods in the inset. 
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FIG. 2. The mean hysteresis-loop area A, defined in Eq. (||), vs. the applied field frequency v 
for T=20K and 100 K. The standard error in the measured mean is smaller than the symbol size. 
The maximum in A, which is the dissipative part of the nonlinear response function, is associated 
with resonance. The temperature-dependent resonance frequencies agree well with those indicated 
in Fig. 1. Histograms of A for several different applied field frequencies at 100 K are shown in the 
inset. 
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FIG. 3. Histogram of the period- averaged magnetization Q, denned in Eq. (||), for applied field 
frequencies of 0.25, 0.14, and 0.1 GHz. Q is a dynamic order parameter for the system. The two 
peaks a the highest frequency indicate a broken symmetry state. 
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